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Abstract 

The ergodic capacity of underlay cognitive (secondary) dual-hop relaying systems is analytically in¬ 
vestigated. Specifically, the amplify-and-forward transmission protocol is considered, while the received 
signals undergo multipath fading and shadowing with non-identical statistics. To efficiently describe 
this composite type of fading, the well-known generalized-iT fading model is used. New analytical 
expressions and quite accurate closed-form approximations regarding the ergodic capacity of the end-to- 
end communication are obtained, in terms of finite sum series of the Meijer’s-G function. The analytical 
results are verified with the aid of computer simulations, while useful insights are revealed. 

Index Terms 

Amplify-and-forward (AF), cognitive systems, ergodic capacity, generalized fading channels, per¬ 
formance analysis. 


I. Introduction 

Underlay cognitive transmission represents one of the most popular spectrum sharing tech¬ 
niques, where secondary (unlicensed) users utilize the spectrum resources of another primary 
(licensed) service. Due to its mode of operation, the transmission power of secondary users 
is limited, such that its interference onto the primary users remains below prescribed tolerable 
levels. However, this dictated constraint dramatically affects the coverage and/or capacity of 
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the secondary communication. Such a condition can be effectively counteracted with the aid of 
relayed transmission. Performance assessment of these systems has been well-investigated in the 
open technical literature to date (e.g., see lH] and references therein). Nevertheless, these works 
assumed non-shadowing environments; a rather infeasible condition. 

In practical wireless communication systems, the signal always experiences composite small- 
scale (multipath) fading and large-scale shadowing simultaneously. The rigorous log-normal 
distribution appropriately describes the latter effect, giving rise to composite models, such as the 
Rayleigh/log-normal distribution. Alternatively, the generalized-iT {Kq) distribution model can 
efficiently describe this composite effect, while preserves mathematical tractability at the same 
time. It is noteworthy that it includes the classical Rayleigh, Nakagami-m and Rayleigh/Gamma 
(i.e., the K distribution) fading models as special types Gl. ll^. 

Yet, only few research works have investigated the performance of cognitive relayed trans¬ 
mission over composite fading/shadowing channels. Specifically, the performance of an underlay 
cognitive relaying system was analytically studied in Q under Kq fading channels, by consider¬ 
ing the decode-and-forward (DF) relaying scheme. The derived expressions therein were provided 
in terms of an infinite series representation. Further, the authors of O studied the scenario of 
underlay cognitive systems with multi-hop/multi-relay transmission under Kq fading channels, 
when the amplify-and-forward (AF) relaying scheme is used. However, the end-to-end (e2e) 
performance was only approximated in that work with the aid of bound expressions, not exact 
ones. Notably, these bounds were tight only in high signal-to-noise ratio (SNR) regions (i.e., 
see [[51 Eq. (6)]). Nonetheless, it should be noticed that the mode of operation used for underlay 
cognitive systems fundamentally supports quite a low transmit power, which contraindicates the 
accuracy of the former bounds in low SNR regions. In addition, identically distributed fading 
channels were assumed in [|51 with common statistics (i.e., equal distance and fading severity 
for all the included signals), which is not always the case in real-life network setups. 

Capitalizing on the aforementioned observations, an underlay cognitive relayed system is 
investigated in current work, where the relay utilizes the cost-effective and computational-efficient 
AF transmission scheme. The signal of each link (from both secondary and primary nodes) is 
subject to independent and non-identical composite multipath fading/shadowing, modeled by the 
Kq distribution. This condition is suitable for most practical applications where the involved 
signals undergo arbitrary link distances with distinct fading statistics. New analytical expressions 
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and quite accurate closed-form approximations for the e2e ergodic capacity are derived, while 
some useful engineering insights are also obtained. The derived expressions are valid in the 
entire SNR region (low-to-high), while they are time-efficient in comparison to other existing 
methods so far (e.g., numerical manifold integrations or Monte-Carlo simulations). 

II. System Model 

Consider an underlay (secondary) dual-hop system where the source (S) communicates with 
the destination (D) via an intermediate relay (R) node. This system operates in the presence of a 
licensed (primary) node (Pr). Moreover, the signal transmission power of the secondary system 
is, in principle, maintained quite low in order not to dramatically affect the reception quality of 
the primary communication in terms of interfering power. To this end, assume that the direct 
communication between S and D is not feasible due to strong propagation attenuation and/or 
severe shadowing, whereas keeping in mind the constrained transmission power regime. Hence, 
the e2e communication is facilitated with the aid of R. 

A. Power Allocation 

The transmitted power of S and R are denoted as Ps = w/ {hspi^f and Pp = w/ \hppj^f, 
where \hspj^f and \hppj^f correspond to the channel gains of S-to-P^ and R-Io-Pr, respectively. 
Also, w represents a power threshold, the so-called interference temperature, which should not 
be exceeded during the transmission of secondary nodes. Such an approach has been widely 
adopted in the open technical literature (e.g., see Cl-llSlI and references therein), mainly because 
it effectively balances performance and complexity. 

In principle, channel state information (CSI) of the links between the primary and secondary 
nodes can be obtained through a feedback channel from the primary service and due to the 
channel reciprocity. CSI can also be captured through a band manager that mediates the exchange 
of information between the primary and secondary networks [|T||. It is noteworthy that when 
w/ \hspf\^ (or w/ \hRpJ^) happens to be higher than the maximal allowable transmitted power, 
say Pmax, power control of the corresponding secondary node may modify w to w' so as 
w’/ \hspj^ = Pmax (or / \hRpJ^ = Pmax) IS Satisfied. This issue is further analyzed in the 


next section. 
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B. Signal Model 

The e2e communication occurs in two consecutive transmission phases, one for each hop. The 
received signal of R at the end of the first phase is given by dr = Iisrx + hr, where i/r, Iirr, x 
and riR represent the received signal, the channel coefficient of the S-to-R link, the transmitted 
signal and the additive white Gaussian noise (AWGN) at R, respectively. Then, this signal is 
amplified with the variable gain G and forwarded to D during the second transmission phase. 
Hence, the overall signal at D is expressed as yo = GliRoyR+nD = GlisRliRox+GliRonR+nD, 
where Iird and no correspond to the channel coefficient of the R-to-D link and the AWGN at D, 
respectively. For notational simplicity and without loss of generality, assume that noise powers 
are identical, i.e., = Nq. 

We retain our focus on CSTassisted AF transmission, since the knowledge of CSI repre¬ 
sents a requisite for the efficient operation of underlay cognitive relaying. Thus, since Pr = 
G‘^{\hsR\^ Ps + No), it can be seen that + ^))-k After some straight¬ 

forward algebra, the e2e SNR reads as 

_ lliRof \hsRf Ps _ 7 i 72 

G^\hRnfNo + No ~ 11 + 12+ G ^ 

where 71 = — and 72 = — “ ' 2 ■ 

1^0 I f^SPji I -^0 I f^RPji I 


III. Performance Analysis 


Ergodic capacity of the entire e2e communication is defined as 


C4iE[log,(l+7.2.)l = fE 


log 2 


(1 + 7 i)(l + 72 ) 
1 + 71 + 72 


— + G2 — G1+2), 


( 2 ) 


where E[.] denotes expectation, Gi = E[log 2 (l + 7i)] with I G {1,2}, G 1+2 = E[log 2 (l+ 71 + 72 )] 
and the factor 1/2 is due to the involvement of two transmission phases. 

Following lemmas will be quite useful for the statistical derivation of Q. 
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Lemma 1: Probability density function (PDF) of 7 ; with I G {1,2} is obtained by 

/ kiimiinjiNo \ ^ ^A-1 

/ X \ kAim^iQiiw I 

^ Tima)Timji)T{ku)Tikji) 

xGg 

where {*,}} = {SR^SPn} for I = 1, while {i,j} = {RD,RPji} for I = 2. Also, 
stands for the Meijer’s-G function |l6l Eq. (9.301)]. Moreover, mu > 0.5 and ku > 0 denote 
the multipath fading and shadowing severity of the ith channel coefficient, respectively, whereas 
higher (lower) values indicate less (more) severe fading channel conditions. In addition, A = 
and 0 = are introduced for notational simplicity, while Vtu = where du 

and au represent the corresponding normalized link distance (with a reference distance equal to 
1 km) and path loss factor. Usually an G {2 — 6}, with an = 2 indicating free-space loss, while 
an > 2 representing suburban to dense urban environments. 

Proof: Since hu is Kq distributed, \hii'^ has a PDF as follows 


kiimuQjiNox 

1 — A — kji, 1 — A — mji 

kjimjiQiiw 

0,-0 


( 3 ) 




kiirriji } 

. J _^ 2,0 

T{mu)Tiku) 


kiimux 

— 


0,-0 


( 4 ) 


For the derivation of dH), the transformation (61 Eq. (9.34.3)] is used into (21 Eq. (2)]. Further, 
it holds that /|;,,^| 2 /|^,^| 2 (a:) = yf^^ufixy) since \hiif and \hjif are mutually 

independent, while /i | 2 (.) is provided in (U) by substituting i with j. Finally, utilizing (7] Eq. 
(2.24.1.1)] into the latter integral and after some straightforward manipulations concludes the 
proof. ■ 


Lemma 2: Moment generating function (MGF) of 7 / yields as 

A 




kiimunjiNo 
kjl PLljl 


T{mu)T{mji)T{kii)T{kji) 


X G 


2,3 

3,2 


kiimukljiNo 

1 — A,1 — A — kji,l — A — mji 

kjlTTlji^ilWS 

0,-0 


( 5 ) 
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Proof: Since = E[exp(— 57 /)], the following integral appears 


^ exp(—sx) 


X G 


2,2 

2,2 


kiima^jiNox 

1 — A — kji, 1 — A — rriji 

kj 1 'kklj 1 ^il W 

0,-0 


dx, 


( 6 ) 


whieh is evaluated in elosed-form with the aid of [|7l Eq. (2.24.3.1)] and ([5]) direetly arises. 
Lemma 3: Partial derivative of is derived by 




xGS 


/ kumuQjiNo \ ‘ 
y kjimjinuw j 


V{mii)V{mji)V{kii)V{kji)s 


A+l 


kiima^jiNo 

—A, 1 — A — kji, 1 — A — rriji 

kjlTTljlQilWS 

0,-0 


Proof: From ®, the following integral appears 


cx 




x^ exp(—sx) 


kiimuQjiNox 

1 

> 

1 

1 

> 

1 

kjlTTljlQilW 

0,-0 


dx. 


Again, using |I71 Eq. (2.24.3.1)], as for the derivation of (l5|), dV]) is extraeted. 
We are now in a position to formulate the e2e ergodie eapaeity. 

Proposition: Ergodie eapaeity of the eonsidered system is expressed as 

A 




kiimunjiNo 

hjlTTljl^^lW 




X g^;4 


kjirrijiniim 


l—A—kji,l—A—mji,—A,l—A 

e-e-A-A 


N 


V=1 


(7) 


( 8 ) 


X (9) 

where Ei(.) stands for the exponential integral funetion Ei Eq. (8.211 .1)],^)^ = —— 

AJM COS ( COS ( j + ■4 j 

and Sy = tan(| cos(^^^^^^^) + j). Ideally, M goes to infinity, but as was shown into |[8l, |0 
(and it will be verified from the subsequent numerieal results), setting AA = 60 provides quite a 
high aeeuraey level. 
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Proof: Based on Ci is obtained by using the transformation of the logarithm funetion 
into the Meijer’s-G funetion [|71 Eq. (8.4.6.5)] and then utilizing ^ Eq. (2.24.1.1)]. Unfortunately, 
this standard teehnique ean not be used for the derivation of G 1 + 2 , sinee the eorresponding PDE 
of 7 i + 72 is not feasible. Nonetheless, based on [jS] Eq. (10)], we are able to bypass the PDE 
involvement and to evaluate G 1+2 as 



( 10 ) 


Still, ([To]) eannot be resolved in elosed-form in terms of standard build-in funetions in well-known 
mathematieal software tools. Thus, using an alternative representation of (fTO]) with respeet to the 
Gauss-Chebyshev quadrature formula (S] Eq. (8)], (fTOl) beeomes 




( 11 ) 


yielding (0). 


Henee, plugging Q and ([7]) into (l9|), a elosed-form approximation of the e2e ergodie eapaeity 
is obtained in terms of finite sum series of the Meijer’s-G funetion. Notiee that G™^”[.] is ineluded 
as standard build-in funetion in most popular mathematieal software paekages; thereby it can be 
calculated quite easily and efficiently for arbitrary fading/shadowing parameters. 

Special case: When transmission power equals P,nax 

Eet 00 , where the kih. hop indicates the asymptotic secondary-to-primary node link 

(i.e., k G {1,2}, k I, conditioned on Vtjk —)■ 0+). The physical meaning of such a scenario 
relies on the fact that when the distance of S-to-P/j or R-to-Pj? link (or both) is considerably 
long, this assumption takes place. Alternatively, this asymptotic behavior can also be modeled 
for relatively shorter distances, but in dense propagation environments (e.g., urban terrestrials 
where aji > 4). Hence, —)■ 00 , reflecting that S> p;^ = Pmax- In other words, the 

conventional (non-cognitive) transmission approach occurs in this case. 







To this end, referring baek to @ and based on ^ Eq. (7)], Ck beeomes 

A 


Ck = 


kikmikNo \ ‘ 

^ik^max. J 


/^4,1 
"^2,4 


" 

” 

>=ik^ik^0 

-A,l-A 

-Pmax 

0,-0,-A,-A 


(12) 


iog(4)r(mjfc)r(fcifc) 

For the more ehallenging C 1+2 parameter in ^ (whieh beeomes in this ease Ci+k or Ck+i), we 
have from [|2l Eq. (4)] that 


-^7fe ('®) 


kikmikNo\^ 2 


exp 


kikiriikNo \ 


^ikPmax^ J \ / 

( kik'kklikP^o\ 

X >V_A+ie 


(13) 


^ikPmax^ ) 

where Wa,/ 3 (.) stands for the Whittaker’s-fE hypergeometrie funetion |j 6 l Eq. (9.220.4)]. Finally, 
following similar lines of reasoning as for the derivation of (|7|), partial derivative of 
yields as 

A-i 




nikPr.j') ( kikTUikNo ^ 

' _GXp I _ I 


2 


,,, f kikmikNo\ 


(14) 


Thus, by appropriately substituting (fT^ . (fT3l) and (fT4l) into Q, ergodie eapaeity is obtained 
whenever 77 -^ > P^ax oeeurs in one or both hops. 

Further, it is interesting to indieate the ease when the latter eondition oeeurs in terms of the 
distanee between the primary and seeondary nodes. By averaging out the presenee of ehannel 
fading based on ^ Eq. (5)], we have that 

' P 

J- m; 


w 


> Pmax ^ d- 


Otjk 


< 


W 


djk > 


1 

^jk 


\h,k? ' ' Pmax “ \ W 

As an illustrative example, let P^ax = OdB, w = —3dB and Ujk = 4 (i.e., a typieal rural 
terrestrial). Then, only when the distanee between the seeondary souree (or the relay) and the 
primary node is greater than 1.19 km, the eorresponding seeondary node(s) may reaeh to the 
maximal aehievable transmission power. 


IV. Numerical Results and Discussion 


In this seetion, the theoretieal results are presented (by setting M = 60) and eompared with 
Monte-Carlo simulations. Table U illustrates an indieative performanee example in terms of the 



















TABLE I 


Number of Terms Required For Convergence LFp To The 4th Decimal Point* 



= d{j,2) 

w/No = OdB 

w/No 

= lOdB w/No = 15dB 

0.05 


49 

51 

51 

0.1 


52 

56 

57 

0.3 


57 

58 

58 

0.8 


58 

59 

60 

= 1, kj 

1 = 4, rrii 

1 = 1, ruj^i = 3, 

a = 4, 

and Pmax/iVo = 20dB. 



Fig. 1. Performance of the normalized e2e ergodic capacity vs. various w/No values. 


included series convergence. Similar behavior is observed for other system parameters. There 
is a good match between all the analytical and the respective simulation results and, hence, the 
accuracy of the proposed approach is verified. In what follows and without loss of generality, 
an identical path-loss factor is used for each link, namely, a = 4. 

In Fig. [H symmetric link distances are assumed between the two hops. Obviously, fading 
severity affects the e2e performance. However, it can be seen that the ergodic capacity is affected 
much more drastically by the link distance of P/? (i.e., dj^i) rather than the fading severity of the 
secondary signal. 

From a different standpoint. Fig. [2] indicates how the ergodic capacity is influenced from 
various (non-symmetric) link distances. Curves in solid lines indicate closer primary-to-source 
distance, while the ones in dashed lines indicate closer primary-to-relay distance. It is clear that 
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Fig. 2. Performance of the normalized e2e ergodic capacity vs. various link distances between the secondary nodes, where 
di^i = d and di _2 = 1 — d. Also, assume common fading parameters for all the involved links, namely, let fc = 2 and m — 1, 
indicating the shadowing and multipath fading severity, respectively. 


when the primary node is nearer (farther) to the secondary source (relay) or vice versa, then 
preserving symmetric distances between the secondary nodes (i.e., when d = 0.5) is not a fruitful 
option. 
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